/**
 * Created by yifangts on 2020/10/15.
 */
const xPI = (3.14159265358979324 * 3000.0) / 180.0
const PI = 3.14159265358979324
const outOfChina = (lng, lat) => {
  if (lng < 72.004 || lng > 137.8347) return true
  if (lat < 0.8293 || lat > 55.8271) return true
  return false
}

const GPS = {
  delta(lng, lat) {
    const a = 6378245.0 //  a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
    const ee = 0.00669342162296594323 //  ee: 椭球的偏心率。
    let dLat = GPS.transformLat(lng - 105.0, lat - 35.0)
    let dLng = GPS.transformLng(lng - 105.0, lat - 35.0)
    const radLat = (lat / 180.0) * PI
    let magic = Math.sin(radLat)
    magic = 1 - ee * magic * magic
    const sqrtMagic = Math.sqrt(magic)
    dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * PI)
    dLng = (dLng * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * PI)
    return [dLng, dLat]
  },
  // WGS-84 to GCJ-02
  gcj_encrypt([wgsLng, wgsLat]) {
    if (outOfChina(wgsLng, wgsLat)) return [wgsLng, wgsLat]
    const [lng, lat] = GPS.delta(wgsLng, wgsLat)
    return [Number(wgsLng) + Number(lng), Number(wgsLat) + Number(lat)]
  },
  // GCJ-02 to WGS-84
  gcj_decrypt([gcjLng, gcjLat]) {
    if (outOfChina(gcjLng, gcjLat)) return { lat: gcjLat, lng: gcjLng }
    const [lng, lat] = GPS.delta(gcjLng, gcjLat)
    return [Number(gcjLng) - Number(lng), Number(gcjLat) - Number(lat)]
  },
  // GCJ-02 to WGS-84 exactly
  gcj_decrypt_exact([gcjLng, gcjLat]) {
    const initDelta = 0.01
    const threshold = 0.000000001
    let dLat = initDelta
    let dLng = initDelta
    let mLat = gcjLat - dLat
    let mLng = gcjLng - dLng
    let pLat = gcjLat + dLat
    let pLng = gcjLng + dLng
    let wgsLat
    let wgsLng
    let i = 0
    const flag = true
    while (flag) {
      wgsLat = (mLat + pLat) / 2
      wgsLng = (mLng + pLng) / 2
      const tmp = GPS.gcj_encrypt(wgsLat, wgsLng)
      dLat = tmp.lat - gcjLat
      dLng = tmp.lng - gcjLng
      if (Math.abs(dLat) < threshold && Math.abs(dLng) < threshold) break

      if (dLat > 0) pLat = wgsLat
      else mLat = wgsLat
      if (dLng > 0) pLng = wgsLng
      else mLng = wgsLng
      i += 1
      if (i > 10000) break
    }
    // console.log(i);
    return [wgsLng, wgsLat]
  },
  // GCJ-02 to BD-09
  bd_encrypt([gcjLng, gcjLat]) {
    const x = gcjLng
    const y = gcjLat
    const z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * xPI)
    const theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * xPI)
    const bdLng = z * Math.cos(theta) + 0.0065
    const bdLat = z * Math.sin(theta) + 0.006
    return [bdLng, bdLat]
  },
  // BD-09 to GCJ-02
  bd_decrypt([bdLng, bdLat]) {
    const x = bdLng - 0.0065
    const y = bdLat - 0.006
    const z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * xPI)
    const theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * xPI)
    const gcjLng = z * Math.cos(theta)
    const gcjLat = z * Math.sin(theta)
    return [gcjLng, gcjLat]
  },
  // WGS-84 to BD-09
  gcj_bd(lnglat) {
    return GPS.bd_encrypt(GPS.gcj_encrypt(lnglat))
  },
  // WGS-84 to Web mercator
  // mercatorLat -> y mercatorLng -> x
  mercator_encrypt([wgsLng, wgsLat]) {
    const x = (wgsLng * 20037508.34) / 180.0
    let y = Math.log(Math.tan(((90.0 + wgsLat) * PI) / 360.0)) / (PI / 180.0)
    y = (y * 20037508.34) / 180.0
    return [x, y]
  },
  // Web mercator to WGS-84
  // mercatorLat -> y mercatorLng -> x
  mercator_decrypt([mercatorLng, mercatorLat]) {
    const x = (mercatorLng / 20037508.34) * 180.0
    let y = (mercatorLat / 20037508.34) * 180.0
    y = (180 / PI) * (2 * Math.atan(Math.exp((y * PI) / 180.0)) - PI / 2)
    return [x, y]
  },
  // two point's distance
  distance(latA, lngA, latB, lngB) {
    const earthR = 6371000.0
    const x =
      Math.cos((latA * PI) / 180.0) *
      Math.cos((latB * PI) / 180.0) *
      Math.cos(((lngA - lngB) * PI) / 180)
    const y = Math.sin((latA * PI) / 180.0) * Math.sin((latB * PI) / 180.0)
    let s = x + y
    if (s > 1) s = 1
    if (s < -1) s = -1
    const alpha = Math.acos(s)
    const distance = alpha * earthR
    return distance
  },
  outOfChina(lng, lat) {
    if (lng < 72.004 || lng > 137.8347) return true
    if (lat < 0.8293 || lat > 55.8271) return true
    return false
  },
  transformLat(x, y) {
    let ret =
      -100.0 +
      2.0 * x +
      3.0 * y +
      0.2 * y * y +
      0.1 * x * y +
      0.2 * Math.sqrt(Math.abs(x))
    ret +=
      ((20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0) /
      3.0
    ret +=
      ((20.0 * Math.sin(y * PI) + 40.0 * Math.sin((y / 3.0) * PI)) * 2.0) / 3.0
    ret +=
      ((160.0 * Math.sin((y / 12.0) * PI) + 320 * Math.sin((y * PI) / 30.0)) *
        2.0) /
      3.0
    return ret
  },
  transformLng(x, y) {
    let ret =
      300.0 +
      x +
      2.0 * y +
      0.1 * x * x +
      0.1 * x * y +
      0.1 * Math.sqrt(Math.abs(x))
    ret +=
      ((20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0) /
      3.0
    ret +=
      ((20.0 * Math.sin(x * PI) + 40.0 * Math.sin((x / 3.0) * PI)) * 2.0) / 3.0
    ret +=
      ((150.0 * Math.sin((x / 12.0) * PI) + 300.0 * Math.sin((x / 30.0) * PI)) *
        2.0) /
      3.0
    return ret
  }
}
export default GPS
